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In  this  paper ,  the  steady,  irrotational ,  subsonic  flow  of  a  gas  around  a 
given  profile  is  studied  in  the  case  of  arbitrary  space  dimension  greater  than 
two.  We  prove  that  the  solution  of  this  problem  exists,  is  unique,  and 
depends  continuously  on  the  incoming  flow.  This  extends  the  previous  results 
of  Bers  and  of  Finn  and  Gilbarg. 
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SIGNIFICANCE  AND  EXPLANATION 


It  is  clear  that  the  study  of  the  gas  flow  around  a  body  moving  at 
subsonic  speed  is  of  primary  importance  in  aerodynamics*  Extensive  studies 
have  been  done  on  the  problem  of  existence  and  uniqueness  of  steady 
irrotational  subsonic  flow  of  a  perfect  gas  past  a  given  profile*  Bers  (1954) 
treats  the  plane  subsonic  flow  using  the  theory  of  quasianalytic  functions. 
Finn  and  Gilbarg  (1957)  deal  with  three  dimensional  flows  with  Mach  number 
less  than  0*7.  In  this  paper  we  prove  the  existence  and  uniqueness  of  three 
dimensional  flow  by  using  and  improving  a  priori  estimates  obtained  earlier  by 
several  authors*  Our  results  can  also  be  extended  to  higher  space 
dimensions#  Furthermore  we  allow  the  flow  to  be  arbitrarily  close  to  sonic 
speed# 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


THREE-DIMENSIONAL  SUBSONIC  FLOWS  AND  THEIR  BOUNDARY 
VALUE  PROBLEMS  EXTENDED  TO  HIGHER  DIMENSIONS 


* 

Guang-Chang  Dong 

The  existence  and  uniqueness  of  steady  irrotational  subsonic  flow  of  a  perfect  gas 
past  a  given  profile  has  been  studied  extensively  in  the  two  dimensional  case.  In  {1], 
Bers  proved  the  existence  and  uniqueness  of  plane  subsonic  flow  around  a  given  profile. 

For  higher  space  dimensions,  few  results  have  been  previously  obtained.  Finn  and  Gilbarg 
[2]  proved  existence  and  uniqueness  in  three  dimensions  provided  the  velocity  was  not  too 
large  (the  maximum  Mach  number  less  than  0.7). 

In  this  paper  we  use  the  idea  of  [2]  together  with  an  improved  a  priori  estimate 
(extending  the  method  of  (4))  to  prove  the  existence  and  uniqueness  of  the  solution  in  the 
three  dimensional  case.  We  also  extend  the  result  to  higher  dimensions. 

In  the  following  we  always  suppose  the  dimension  of  space  is  n(n  >  ?),  and  use  the 
summation  convention 


a  x.  *  T  a  x  , 
ii  .  .  li 

1*1 

and  denote  the  vector  (x1,...,xn)  by  x. 

The  steady  irrotational  gas  flow  in  n  dimensional  space  can  be  described  by  the 

velocity  potential  v(x),  satisfying  the  equation 

jp-  (p  -  0  (•> 

dx  3  X 

i  i 

where  p  represents  the  density  of  gas,  which  is  a  given  positive  function  of  velocity 
q,  where 

q 

( 1 )  can  be  written  as 


L  2  3v 

/Eu  . ,  u  .  *  tr — 
l  l  3  x 


3  v 

3i j  3x . 3 x . 
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a.  .  *  p<$.  .  +  2—  u  u  ,  5  . 

i}  q  i  j  ij 


y  (i-j) 
0  U+j) 


(3) 


Because  (1)  and  (2)  are  rotational  invariant,  the  eigenvalues  X . .  of  the  quadratic 

i  n 

form 

can  be  obtained  by  letting  u1  =  q,  u2  =...*=  un  *  0,  i.e. 


X,  =  p  +  p'q.  \2 


Xn  =  p 


So  that  if 

P  +  p'q  >  0  (4) 

at  this  point  the  quasilinear  equation  is  of  elliptic  type,  and  the  flow  (1)  is  subsonic  at 

this  point.  If  in  some  region  ft  we  have 

inf ( p  +  p'q)  >  0  , 

ft 

then  the  flow  (1)  is  a  subsonic  flow  in  ft. 

Assume  that  the  function  p  =  p(q)  is  a  positive,  three  times  continuous 

differentiable  functions,  defined  in  the  interval 

0  <  q  <  q  „ .  4  ®  , 

£im 

and  assume  that  a  number  a  exists  (we  called  it  the  critical  velocity)  0  <  q  <  q  .  , 

c  £im 

such  that 


p  +  oo 


da 


(  pci )  >  0  ( 0  <  q  <  q  ) 


(5) 


For  definiteness,  assume  that  when  q  >  ,  the  left  hand  side  of  (4)  is  non-positive. 

Hence  (2)  is  subsonic  if  and  only  if  0  <  q  <  .  We  also  assume  that 

p'(0>  =  0  .  (6) 

In  gas  dynamics,  the  density  of  a  perfect  gas  is  given  by  the  dimensionless  formula 


o'  (l  -  2  <1* 


1 

y-1 


(1  <  y  <  2,  y  is  a  constant) 

Irr 

fi  -  x 

2 


y-1 


,  , .  Y-1  2  N ’  '  / ,  Y+1  2 

p  +  p'q  =  i  -  hr  q  J  l1  “  hr  q 


2 


hence 


/-£- 

/y+t 


Him 


/-i- 

Y-1 


The  flow  Is  called  subsonic  when  0  ^  q  <  q  ,  which  agrees  with  the  discussion  above. 

c 

The  profile  T  is  a  bounded  closed  surface  in  n  dimensional  space  (it  can  be  some 
closed  surface  also) .  Assume  that  a  constant  exists  such  that  0  <  tq  <  1  and 


(2  +  T  ) 

r  e  c  (7) 

The  region  outside  T  (not  containing  the  point  at  infinity)  is  denoted  by  ft.  We  only 
consider  those  satisfying  the  condition:  any  closed  curve  in  ft  can  be  deformed  into  a 
point  without  touching  T • 

The  flow  around  T  means  that  the  solution  of  ( 1 )  in  ft  satisfies  two  conditions  as 
follows-  The  boundary  condition 


Hi 


0  (N  is  the  interior  normal  of  ft  on  D 


and  the  condition  of  uniform  flow  at  ft  (called  uniform  incoming  flow) 

00  00 

grad  i j  x _  =  u  (u  is  a  constant  vector)  • 


(8) 


(9) 


Usually  for  the  problem  of  flow  around  a  given  profile  we  give  the  magnitude  and 
direction  of  incoming  flow.  Without  loss  of  generality  we  can  assume  that  the  direction  of 
incoming  fiow  is  the  positive  x^  axis  (otherwise  after  an  axis  rotation  we  can  reach  this 
a ituation) ,  i. e. 


u°°  =  (U ,  0 , .  •  • ,  0  )  the  const  U  >  0  is  given  •  (10) 

A  slightly  different  problem  is:  Given  the  direction  of  incoming  flow  but  not  the 
magnitude,  i.e. 

u  *  (U,  0, • -  - ,0 )  U  >  0  is  not  given  (11) 

and  given  the  maximum  value  of  speed  instead,  i.e. 

sup  q  “  Q  is  given  .  (12) 

ft 

We  shall  study  the  existence  and  uniqueness,  find  the  flow  and  other  properties  of  the 
solution  of  the  above  two  kinds  of  flow  problem  in  the  subsonic  range,  i.e.  when 

Q  <  %  •  (13> 

$1.  Preliminary  study  of  the  ’  inoar  problem 
Consider  the  following  problem 


-3- 


(1.1) 


Lm  -  Vx)  aSr +  bitx)  fr  ■ f(x)  (x  e  n'  bii  -  v 


3t 


3ulr 


t^  (N  is  the  interior  normal  of  ft  on  f) 


(1.2) 


t|  *  0 


(1.3) 


where  the  given  functions  b^ f  b^,  f  c  C(0)(ft),  e  C(0)(D,  and  a  positive  constant  o 


exists  such  that 


r  2  .  .  1  _  2 
Ola.  <  b.  a  a.  <  -  Da 
i  i]  i  ]  a  i 


(1.4) 

Further  conditions  on  b^,  f  are  given  below.  Assuming  that  the  solution  of  problems 

(1.1)#  (1.2),  (1.3)  exists,  and  we  deduce  some  a-priori  estimates. 


Assume  that  a  fixed  small  positive  number  x  such  that  x  <  rQ  (the  definition  of 


see  (7))  exists  such  that 


lb. .(x)  -  b  (•) |  4  K  r  T(x  e  ft) 

1  3  ^  J  T  X 

lb  (x)  -  b  ( y )  |  4  K  r"^|x-y|T  (x,  y  e  ft) 
aj  xj  *  xy 


'V  <  Vx 


1  —  T 


(1.5) 

(1.6) 

(1.7) 


where  the  constants  K^#  K2,  ...  depend  on  o,  x,  F  only.  Moreover,  assume  f  is  a 
bounded  function  satisfying 


| f |  <  K  max  | f | r  n  T  (1.8) 

4  a  x 


where 


r  *  d(x,  D,  r  -  min{r  #  r  ) 
x  xy  x  y 

d(x,  T)  is  the  distance  between  x  and  T •  We  have 
Lemma  1 . 

It |  <  K  (max  |t  |  +  max  |f|J  . 

5  r  0  ft 

Proof.  Let 


(1.9) 


(1.10) 


t  *  max  1 1  |  ♦  max  |  f  | 

r  ft 


(i.ii ) 
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The  uniqueness  of  the  solution  of  (1.1),  (1.2),  (1.3)  is  easily  deduced  from  the 


maximum  principle  [5].  Thus  for  7  *  0  we  have  7=0,  and  in  this  case  (1.10)  is 


When  7^0,  let 


*  V“’  aJpsq 


(1 


R  =  /bTT ( *)x  .  x7  . 
i  j  13 

B  («)  is  the  algebraic  complement  in  the  matrix  (b. .(«)).  We  have 

ij  13 


<  1 


2-n-  ~  -n-  ^ 

L  (R  )  =  j  (n-2+  1)R  ‘ 


(  1 


When  R  >  K  ,  where  KA  is  a  sufficiently  large  constant,  from  (1.5),  (1.7),  (1 
o  ^ 

(1.14)  we  have 


2-n-  —  -n-  — 

l(r  2)  -  ^  ~  (n-2+  ~)R  2  +  0(R~n'T)  >  0 


7 


(  1 


Without  loss  of  generality  we  can  assume  that  the  origin  lies  inside  F,  then 

min  R  >  0.  Take  a  ~  a(a#T#D  comparatively  large,  by  using  (1.5),  (1.7),  (1.13)  we 

xeft 


L(R~a)  *  a{a+2)R’0r4b.  .B.  ,(»)x  3.  (®)x  -  aR”0^2*?.  .  ft .  .  ( »)  - 

lj  x]  k  ]h  n  i] 


-  OR  Vik(”)Xk  >  K1R 


-a-2 


(  1 


From  (1.15),  (1.16)  we  can  select  a  constant  such  that 


2-n-  2  -  f 

L(R  +  KRR  a)  - -  >  0  (x  e  2) 

7 


(1 


/  ^  =  y  ;  y°(R  2  +  krr  n  . 


From  (1.1),  (1.17),  (1.1ft)  we  have 


(  1 


L(f  )  <  0,  L(72)  >  0 


(  1 
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true. 

.12) 

.13) 

.  14) 

.8), 

.15) 

have 

.16) 


.17) 

.18) 

.  19) 


I 


\ 


From  (1.3)  we  have  ^t®)  15  0,  hence  min  f  <  0,  this  minimum  can  only  be  taken  at 

a 

»  or  on  r  by  (1.19).  Similarly  max  4^  >  Or  and  this  maximum  can  be  taken  at  ®  or 

on  F.  Hence  we  have  four  cases  as  follows. 

Case  1.  min  4^  =  ^(Pq)  <  0,  PQ  e  ?t  max  *  4'2(Qq)  y  °*  c  v  • 

( 2 ) 

Extend  the  method  of  (6)  to  prove  that  (1.10)  is  true.  Because  T  e  C  ,  there 

exists  a  positive  constant  y  =  y(P),  such  that  for  any  point  P  e  T  we  can  draw  an 

exterior  tangent  sphere  with  radius  y  which  lies  entirely  inside  ft  except  point  P. 

Draw  the  exterior  tangent  sphere  with  radius  y  at  PQ,  denote  its  center  by 

Consider  the  function  _ 

-hP  P2  2 

V  (P)  »  f^P)  -  min  V  -  kf  (e  -e  7  )  (P  e  ft)  (1.20) 

ft 


where  k#  h  are  positive  constants  we  may  choose  at  our  disposal.  Take  h  *  h(a)  large 
enough,  from  (1.17),  (1-20)  we  have 

L(Y3>  <0,  ^  <  Fjp  <  y  .  (1.21) 

From  (1.20)  we  have 


If 


4^  >  0,  when  P^P  =  y  and  H^fPg)  =  0 


(1.22) 


<P  >  0  when  P^P  *  *  (1.23) 
is  true,  then  from  (1.21),  (1.22),  (1.23)  we  see  that  min  4^  <  0  when  ^  <  P^P  <  y  can 
not  he  true,  hence  min  =  From  (1.2),  (1.1B),  (1.20)  we  have 


3 

3N 


U  *  Vpo’  ' 

0 


2-n- 


v'd,| 


(1.24) 


Pick  the  constant  kQ  so  that 


-h  y 

2yK0he  T  >  1 
then  (  1.24)  can  not  be  true  when  we  take 


max|i5(R  2"KflR'c‘>l 

K  ~  Kq ,  hence  (1.23)  can  not  be  true,  in  other 


words,  we  have  P2  satisfies  p<jP2  °  2  suc^  that  Y^lP^)  <  prom  (1*20)  we  have 

(1.25) 

<£  K.  \J  £. 

Similarly  there  exists  a  Q2  satisfying 

V  1v 

(1.26) 

Because  of  (1.6)  we  can  apply  Barnacle' a  inequality  [7]  to  the  nonnegative  function 


VV  "  min  yi  4  K,(/'  2  <  dtV  V  4  ¥  • 

ft 


max  12  -  V2 (Q2  >  <  K^0,  j  <  d(C2'20)  <  F  * 


-  min  Y^  and  obtain 


Y  AQ  )  -  mm  Y  <  K  [Y  (P  )  -  min  Y  +  max|L(Y  )|] 
^  ^  jj  Tai  1  2  g  1  ft  1 


(1.27) 


Combine  (1.11),  (1.18),  (1.25),  (1.27)  we  have 


W  -  mf  *1  «  *13* 


(1.28) 


From  (1.18),  (1.26),  (1.28)  we  have 

ma  Y  -  min  Y  <  max  Y 


min  Y  <  (K  +  K  )Y 
*  ft  11  J 


ft  ft  ft 

From  (1.3)  we  have  max  Y  >  0,  min  Y  <  0,  so  combining  the  above  expression  and  (1.11)  we 

ft  ft 

obtain  (1.10). 


Case  2.  min  Y  =*  Y  («)  =  0,  max  Y 

ft  1  ft 


VV  >  0,  C0  e  r 


By  lemma  2,  there  exists  positive  constants  R^(0,t,r)  and  K14,  such  that  the 

Harnack  inequality  for  the  positive  function  R2  ^Y,  (the  definition  of  R  see  (1.13)) 


for  any  >  RQ 


holds: 


msx  (R2  nY^  )  <  K14  min  ( R2  nY^  )  +  max  |  f ) 

R=RQ  R=R,  « 


(1.29) 


Proof  of  (1.29)  see  lemma  2.  From  (1.29)  we  have 

max  T  -  r"“2  Max  (R2-n't'l  )  <  K)4r"_2  max  (R2_n 
r=rq  R=R0  R"R, 


V  ♦  k16r0"-V  - 


=  K  [  — ' 
-  1  ^  c  J 


From  (1.19)  and  Y 


^  1  |w 

0,  combining  1 

max  Y 

RsR 

*  max  Y  <  ma 

R=R1  R»i 

1.31)  we 

have 

Rn  n*2  o  n 

i4vR;)  max  *i  +  k,5ro’  * 

1  R=R  ^ 


16 


(1.30) 


(1.31 ) 
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11  -  Ku^n’2]  -  ♦  «CV° 


/  RCK  1 

Take  R.  such  that  K.  I  —  I  =  —  we  have 
1  14  R,  2 


max  V1  <  X  Y 
R=R0 


0 


n .  32 ; 


(1,32)  is  similar  to  (1.25).  We  can  get  (1.26),  (1.27),  (1.28),  (1.10)  by  the  similar 
process  as  in  case  1. 

Case  3.  min  Y  =*  Y  (P  )  <  0(P  t  D  ,  max  Y  =  Y  <«)  =  0. 

5  S 

Obtaining  (1.10)  is  similar  to  case  2. 

Case  4.  mm  Y  »  0  =  Y  (°°),  max  Y  =  0  =  Y  (<*»)  . 

5  1  1  5  2 

Prom  (1.18)  we  have  Y  =0,  contrary  to  the  hypothesis. 

Lemma  1  is  thus  proved* 

Lemma  2.  Let 

Hf  >  0  ( x  £  n)  (1.33) 

Prove  that  (1.29)  is  true. 


Proof.  Assume  that  the  region  R  >  R  (the  definition  of  R  st. e  (1.13))  lies  m 

1  O 


S?*  Without  loss  of  generality  we  can  assume 

t  1 . 34  > 
M  .35) 

Otherwise  after  a  linear  change  of  independent  variables  and  multiplying  the  unknown 
function  by  a  suitable  constant  we  get  (1.34)  and  (1.35).  Hence  R  -  |x|.  Apply  the 
inversion  transformation 


max  J  f |  *  1 
ft 

b .  (*»)  *  <5 
A3  13 


and  let 


(1.1)  becomes 


(i  =  1 , 2 , .  • . ,  n ) 


I  xf 


2-n- 


Yi  -  [ R  -  R 


2v 


(1.36) 


(1.37) 


Uv)  =  BijVXiXj  +  Bivx1 


Bv  =  F  (  I X  |  <  — }  . 

*"18 


Notice  that  the  matrix  is  not  the  same  as  in  (1.23).  From  (1.4)  and  (1.36)  we  see 
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L ( 1 )  =  B  <  0,  L  {  1  +  2 | X | 2  >  >  0  .  (1.46) 

We  need  to  determine  continuous  functions  K+(X,  Y),  K_(X,  Y)  in  the  range 
|X|  s  (y|  ~  (X+Y)  satisfying  the  following  three  relations: 

i.  UK+)  >  0,  L(K_)  «.  0  (0  <  IXj  <  K^q)  (1.47) 

n.  lim  /  K  (X,  Y)g(Y)dS  =  g(Y  )  (1.48) 

X+Y  1  ° 

0 

where  g(X)  is  any  continuous  function  and  Y q  is  any  fixed  point  satisfying  |Yq|  =  K->q. 


< 


iii.  There  exists  constants  1 #  *22  9uch  th*t 

K 

K+(X,  Y)  >  K21.  K_(X.  Y)  4  K22  (|X|  <-f-.  |T|  -  K2q)  .  (1.49) 

If  K+  ( X ,  Y ) ,  K_(X,  Y)  exist ,  then  combine  (1.44),  (1.45),  (1.47),  (1.48),  (1.49), 
and  apply  the  method  of  (5],  (c.f.  in  $4  of  [5])  which  extends  theorem  1  to  the 
inhomogeneous  equations.  This  is  the  result  that  we  need.  Hence,  positive  constants  K23' 
*24'  *25  ex*st  SLlch  that 

K 

K23v(0)  -  K24  4  v(X)  4  k"3  v(0)  +  K25  ( | X |  <  ~y)  ,  (1.50) 

(1.29)  follows  from  (1.34),  (1.37),  (1.50). 

We  can  construct  K+  and  K_  by  only  altering  the  method  of  [ 5 ]  a  little  as  follows. 


H  =  (K2  -X2) fBf . (Y) (X.  -  Y  ) ( X .  -  Y.)]  2 

20  1]  i  1  3  3 

where  B*^  are  the  algebraic  complements  of  .  Take  K_  *  f_(H),  then 

L(K_)  =  B  (X)Hx  Hx  (f”  -f  AfM  ♦  Bf_ 


A  =  [B.  -(X)HY  v  +  B,Hy  )  /  (B .  ,  (  X  )HV  Hy  ) 


i  j  XiX 


i j' A,nX .  X  i  * 


Apply  the  estimate  method  of  [5]  (the  method  of  obtaining  (30)  in  it)  and  (1.6), 


(1.39)  and 


I-' 

B.  »  0  ( J  X|  )  =  0  (  | X-Y  |  -  K 


deduced  from  (1.40),  we  have 

i  -  i 

I  A  |  4  K2(.<  |X-Y|n_1  +  T  +  I  X-Y  |  "  |  X-Y  I  -  K20  |  2  )  . 

And  then  determining  f  etc.  (these  steps  are  similar  to  [5]),  at  last  we  get  K_(X,  Y) 
satisfying  (1.47),  (1.4R),  (1.49). 

1 

2 

Through  the  t ransforma t ion  K+  *  ( 1  +  2  I  X I  ) K* ,  (1.47)  becomes 

L*(K*)  >0 


Applying  (  1.46)  we  have 


L*  ( v )  -  B  v  +  B*  v  +  B*  v 

i  ■)  X  .X  ,  l  X  . 

11  i 


B*  =  LM  +  2(X|2)/(1  +  2  f  X  |  2 )  >  0  . 
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Let  K+  *  f  +  { H  )  •  the  selection  of  f+  is  similar  to  [5].  The  proof  of  the  existence 
of  K*  such  that  K*  satisfying  (1.47),  (1.48),  (1.49)  is  similar  to  K_.  Hence  the 
proof  of  lemma  2  is  finished. 

Lemma  3 •  We  have 

!¥ I  <  K  mind,  r2  n)(max|f  |  +  max|f|)  (1.51) 

2  r  °  5 

under  the  conditions  (1.4),  (1.5),  (1.6),  (1.7),  (1.8). 

Proof.  Because  of  lemma  1,  we  only  need  to  prove  the  validity  of  (1.51)  when  r 


large  enough.  From  (1.13),  (1.14)  when  R  >  (^g  lar<?e  enough)  we  have 


L(R2~n  -  R 


-  7 

.  JL 


k5y 


T 

-n-  - 

-  -  -  n  -  2  +  ^)R  +  0  (R“n“T)  <  0  . 


.  2-n-  ~  „  2-n-  ^  . 

(R2”n  -  R  2)/(K2en  “  K28  2)^  ±  Y 


Hence  the  function 

„  2-n-  ~r  ^  2-n-  i  . 

^  -n  J  J  —  0  fl 

(1.52) 

can  not  take  on  a  negative  minimum  value  in  ^g  <  R  <  ®.  From  (1.3),  (1.10)  we  see  that 

when  R  *  oo  and  R  ■  K28  < 1  .52)  is  nonnegative,  hence  (1.52)  is  nonnegative  in 

K  „  <  R  <  oo.  This  proves  the  lemma. 

28 

(  k )  (  k) 

Denote  the  k'th  order  derivatives  of  f  by  D  Vf*  When  f  e  C  (fi)  or 
C{k+T)(Q),  let 

Mm  k(T)  =  lub  rm+k|D(k>Y(x)  |  (1.53) 

m'  xeQ  X 


m  d)  a  lub  r  |D  ?(x)  -  D  y(y ) 1/ |x-y | T 

m,k  +  x  _  xy 

x,vefl 


(1.54) 


where  m  is  a  nonneaative  integral.  For  the  meaning  of  symbols  r^,  rx^  see  (1.9).  The 
symbol  of  lub  in  (1.53),  (1.54)  means  first  to  take  upper  bound  for  all  k's  order 
derivatives.  Let 


m ,  k 


y 

1=o 


M  .  ( T) ,  I'M  .  ^  =  I'M  +  M  ( T ) 

m,k  m,k+x  m,k  m,k  +  T 


-1  1- 


When  Y(x)  e  C(k)(0)  and  IH 


is  finite,  denote  by  ¥(x)  e  C 


Similarly  we  define 


Lemma  4,  Assume  that  (1,4),  (1*6),  (1*7)  and 


IVX>  "  bi(y)l  <  K29|x_y|Trxy~T 


are  valid,  and  assume  that 


Let  be  the  solution  of  (1.1),  (1.2),  (1.3)  and  ¥  c  C  ,  then  we  have 

n-2 , 0 

V  e  C  _  _ .  ,  and 
n-2 , 2+T 

I y  1  ,  ,.T  <  K__[M  -(f)  +  |f|  J  .  (1.57) 

n-2,2+T  30  n-2,0  n,t 

Proof.  Prom  the  interior  estimate  [9]  we  obtain  (1.57). 

-  -  ( k ) 

Consider  a  bounded  region  ft^  satisfying  ft  o  ftQ  T .  when  V  c  C  (ftQ  )  or 

C  C(k+T)(Sl0),  let 

M  (*)  -  lub  |D<k)»(K)|,  M  ...IV)  =  1  ub  |D(k)f(x)  -  D(k)5'(y)  )/|x-yf  , 

xtn0  k  x,y££l0 

|ylk  =  !  M  <*>,  I'M  =  mk  *  Vim  * 

3*0 

When  f(x)  e  C*** (ftQ )  and  |H,|^  is  finite,  denote  by  ¥(x)  c  C^,  similarly  we  define 


Concernng  the  function  defined  on  T,  let 

n  II  -  lub  lub|D(h)('i'  )|,  ll¥  I  *  If  U  +  lub  |D(k,4'n(x) 

0  k  ...  r  00  k+T  Ok  r  0 

h«ik  xer  x,y cT 

-  D(k)4‘  (y)  i/|x~y|  X 

k  { k )  U 

where  in  lub  the  D  ¥g(x)  and  D  4^(y)  restrict  to  the  derivatives  of  the  same 

parameter. 

Lemma  5.  Apart  from  the  hypothesis  of  lemma  4,  let  b^,  b^,  f,  VQ  satisfy  the  following 
conditions 


^bij(x)  ’  bij(y) ^  *  K31lx~yl  *x'  y  C  V 

|b.(x>  -  b^ ( y ) |  <  K32lx~y|T 


And  let  ¥  e  C 
and 


(  2  +  t) 


fQ  e  C(1  +  T)  (D  . 

(Q)  be  the  solution  of  (1.1),  (1.2),  (1.3),  then  we  have  7  e  C 


(1.61) 


Proof.  See  [10] 
Let 


4  *33IVf>  +  lflT  +  #Vl  +  T] 


M  (7)  =»  lub  max(  1 ,  r™** ) | D*k * 7( x) | 

m,  K  _  x 

xen 


2+t' 


(1.62) 


M 


m,k+T 


lub  mil,  rm+k+T)|D(k,'Kx)  -  D(k)  y(  y)  1/ 1  x-y  |  T 


X,  Y  efl 


llfl  ,  ■  y  M  AV),  I7H  _ 
m,k  m,j  m,k+x 


178 


m,k 


M  (7) 

m,  k+i 


Summing  lemmas  3,  4,  5  up,  we  have 

Lemma  6.  Under  the  conditions  (1.4),  (1.5),  (1.6),  (1.7),  (1.8),  (1.55),  (1.56),  (1.58), 

(1.59),  (  1,60),  (  1.61  )  let  7  £  Cf^(ft)  be  the  solution  of  (1.1),  (1.2),  (1.3),  we  have 

178  „  ,  <  K.JIfl  +  )  .  (1.63) 

n- 2 , 2  +  t  34  n,x  0  1+t 

Theorem  1.  If  the  conditions  of  lemma  6  are  satisfied  then  the  solution 
7  £  C*°*(ft)  n  C(2)  (0)  of  (1.1),  (1.2),  (1.3)  exists  and  is  unique.  7  e  C(2  +  T*(r>),  and 
satisfies  the  estimate  (1.63). 

Proof.  The  uniqueness  of  the  solution  follows  directly  from  maximum  principle.  The 
existence  of  the  solution  can  be  obtained  by  applying  the  continuity  method,  i.e.  to  solve 
the  equation  with  parameter 

[M-9)L  +  8L]  7  =  9f  (0  <  0  <  1)  (1.64) 

with  boundary  conditions  (1.2),  (1.3),  where  the  definition  of  operator  is  defined  by 

(1.12).  when  9=0,  through  a  linear  transformation,  the  equation  (1.64),  (1.2),  (1.3) 
becomes  the  Laplace  equation  in  an  exterior  domain  with  oblique  derivative  given.  The 
existence  of  a  solution  7  and  7e  C^2+T^(P.)  can  be  obtained  bv  the  method  of  integral 
equations,  see  [5],  Applying  lemma  6  to  extend  the  solution  by  increasinq  the  parameter 
0  tint  i  1  9  :=  1 . 

The  theorem  is  thus  proved. 
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§2.  The  existence  of  solution  around  a  given  profile. 

Theorem  2.  The  solution  of  problem  (1),  (9),  (9),  (11),  (12),  (13)  exists. 

Proof.  First  we  alter  the  function  p  a  little  as  follows:  Take  Q 
for  the  definition  of  qc,  Q  see  (5),  (12).  Let 

s  =  s(q)  =  j{1+th[(Q-Q)(2q-Q-Q)]/[2(q-Q)(Q-q)]} 

(  ») 

then  s ( q )  e  C  when  Q  <  q  <  Q  and  increases  monotonously*  and 

s(o)  =  S'(Q)  *  sM(Q)  s’(Q)  *  s  '*  ( Q )  =  ...“  0,s(Q)  =  1 


{p  (0  <  q  <  Q) 

p( 1-s)  +  es  (Q  <  q  <  Q) 

e  (Q  <  q  <  ») 


where  the  constant  e  *  max  p,  it  is  easy  to  see  that  the  smoothness  of  p  is  the  same 
(3)  - 

as  p,  i.e.  p  e  C  ,  and  p  >  0.  Applying  (5)  we  have 

p  +  p *q  *  ( p+p'q) ( 1-s)  +  es  +  q(e-p)s*  >0  (Q  <  q  <  Q)  (2.2) 

Consider  the  equation 

jtrp  If.'  =  0 

i  l 

or 

a.  .(Dv5)  -3  -f  =  0  (2.3) 

11  3x. 3x . 

i  3 

where 

—  —  n ' 

a  .  .  ( Ds? )  =  p6.  .  +  u.  u  . 

H  13  q  i  3 

The  eigenvalues  of  its  correspondent  quadratic  form  a^(rv)a^cr  are 

A1  =  p  +  p*q,  x2  ** •**  Xn  =  p  * 

From  (2.2)  we  see  that  all  eigenvalues  are  positive,  i.e.  when  0  *  q  <  «  (2.3)  is  always 


an  elliptic  equation,  and  a  positive  constant  o  ~  o(Q)  exists  such  that 

2  12 
gla.  <  a  .  (Dtf)ct.  n ,  <  —  In. 
l  i3  i  j  o  i 

Next  we  prove  that  the  solution  of  (2.3),  (9),  (9),  (11),  (12)*  (13)  exists.  From 


(13)  and  (2.1)  we  see  that  it  also  is  the  solution  of  (1),  (8),  (9),  (11),  (12),  (13). 


Hence  concerning  the  proof  of  theorem  2,  without  loss  of  generality  we  can  assume  the 
quasilinear  equation  (2)  is  uniformly  elliptic  when  0  <  q  <  ®,  i.e.  there  has  positive 
constant  o  *  a(£)  such  that 

at  or  <  a  (D  ^j)a  a  <  “  Za2.  .  (2.4) 

i  ij  i  j  o  i 

Let  the  function  vp(x)  satisfy 

*  e  C(1+T)  (fl+D  (2.5) 

D^(®)  -  (U,  0,...,0) 

U>(*)  -  UxtJ  -  0  (2.6) 

|D^Kx)  -  D^j(®)  |  <  c-jrxT 

|D^(x)  -  D^?(  y)  |  <  C  min(  1 ,  r  T)|x-y|T 
t  xy 

where 


T  ■  t(Q,D  <  (the  definition  of  tq  see  (7))  (2.7) 

is  a  positive  constant  which  shall  be  determined  later,  U  is  a  positive  constant  not  to 
be  fixed,  (for  the  meaning  of  r x,  rxy  see  (1.9)),  C j ,  C2  are  positive  constants.  The 
norm  of  ^(x)  is  defined  by 

I “  |U|  +  max/^~Dx)  t  inf  C  +  inf  C  . 

s 

It  is  easy  to  see  all  functions  $  form  a  Banach  space  E. 

We  wish  to  find  the  solution  of  the  following  equations: 


32Y 


(2.8) 


aSlr  ’  ‘cos(N<  Vr 


(2.9) 


f|  =  0  . 


(2.10 ) 


Let 


)  (  x )  *  eQlf(x)  +  x  t  ]  /max  |  grad  [  ¥  ( x )  +  x }]  |  (0  <,  0  <  1) 


(2.11  ) 

*  a  1 

where  the  definition  of  g  see  (12). 

By  theorem  1  we  have  0  e  E,  hence  a  functional  in  E  is  defined  as  follows 

$  =  T(*,  0)  (0  *  0  <  1)  .  (2.12) 
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Now  use  Leray-Schauder  degree  theory  (11]  to  prove  that  for  any  0  <  9  <  1  the 
solution  of  functional  eguation 

<p- T(^j  ,  9)  -  0  (2.13) 

exists.  By  theorem  1  and  (2.8),  (2.9),  (2.10),  (2.11)  we  see  that  when  9  =  1  the 
solution  of  (2.12)  is  a  solution  of  (1),  (8),  (9),  (11),  (12),  (13). 

Let  us  verify  the  three  conditions  of  Leray-Schauder  degree  theory, 
i.  From  theorem  1  we  have 


HI  ,  0jL  <  K  I**! 
n-2,2+T  34 


(2.14) 


Hence  T  is  completely  continuous.  T  is  also  closed,  i.e.  when  <0  e  E,  *  t  E  and 
II  ^  -  *1+0  ( k  -►  «»)  the  corresponding  solutions  satisfying  B II  >0.  If  this  is 
not  true,  i.e.  there  is  a  subsequence  {k* }  and  positive  constant  5  such  that 


-  *1  >  6  . 

k' 


(2.15) 


When  kf  is  large,  from  (2.14)  we  have 


H.  .  »  .  _  <  *..<1  ♦  1*1)  ' 

k'  n-2,2+t  34 


hence  {Hf.  ,  }  is  compact  in  II . . .  I 

n-2,2  +  j 

(kH)  of  {k'}  such  that 


i.e.  there  exist  and  a  subsequence 


32Vt 


From  a  .  .  (D*  , ) 


i  i  k"  3x  .  3x  . 

i  3 


Hh  -  fJ 

k  n-2,2  ♦  l 


-  0  taking  limits  we  have 


0  . 


^Ti7-=  0  • 

i  j 

Similarly  we  have 

3  y 

W\r  *  ‘cos(N'  Vr'  M-  =  0  • 

Hence  from  the  uniqueness  of  the  solution  of  (2.8),  (2.9),  (2.10)  we  have  -  H1,  hence 


!♦,  „  -  #|  +  0  . 
k 

This  contradicts  to  (2.1^),  therefore  T  is  closed. 

From  (2.11)  we  see  that  T  is  continuous  uniformly  in  9,  hence  the  degree  of  (2.13) 
is  indenendent  of 

It.  When  8  s  o,  <t>  *  0 ,  i.e.  (2.13)  has  only  the  solution  *  =  0.  Therefore,  the 


degree  of  (2.11)  is  equa1  to  1. 
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iii.  To  prove  the  solution  of  (2*13)  is  bounded  in  E,  in  other  words,  to  prove  that 
the  solution  of  (1),  (8),  (9),  (11)  and 

max  q  *  0Q  (0  <  6  <  1 )  (2.16) 

5 

is  bounded  in  E,  which  shall  be  proved  in  the  following  two  sections,  first  notice  that 
from  (2.16)  we  have 

max  q  <  Q  (2.17) 

5 


§3.  The  interior  and  boundary  estimates  of  solution. 

Let  v?  e  E  satisfy  (1),  (8),  (9),  (11),  (2.17).  Consider  (2)  as  a  linear  equation  in 
s ot  from  (2.5)  and  the  Schauder  estimate  (notice  that  we  have  (7)  and  (2.7))  we  have 

■f  e  c(3+T)  <fi )  n  c(2+T)  <n+n  . 


3  3  u  3s? 

r —  (a.  t — )  =0  (u  is  the  abbreviation  of  t —  =  u, 
dx  .  l  j  dx  .  dx,  n 

i  3  h 


where  the  definition  of  see  (3). 


(3.1) 


Assume  the  sphere  | x-x^ I  *  V  lies  in  Q  and  0  <  0  <  1.  Let 


1  lx-*0l  <  PCI-8) 

C(x)  =  £(x,  li,  0)  =  \  ^  (U  -  |x'x0l)  n(1-0)  <  I x- xQ |  *  u 

0  |x-xQ|  >  y  * 


(3.2) 


Multiply  (3.1)  by  £  (x)u  and  integrate,  we  have 


fc  a, .u  u  dx  *  -2/a. .uu  CC  dx 

4  n  v  y  ;  1-1  v  v 


Xi  X. 


From  this  formula  and  (2.4)  we  have 
2  .  ,  .  2 


il  x .  x  . 
1  3 


oj^igrad  u)^dx  <  L^/teC^grad  u)2  +  —  u2(grad  Q)Z] 


dx 


In  this  and  the  following  sections  the  constants  L ^ ,  L 2, 


depend  on  Q,  V  only.  c  is 


an  arbitrary  positive  constant.  Take  eL^  =  —  we  have 

/^^(grad  u)2dx  <  L^/u^Cgrad  C)2dx 


(3.3) 


Take  0  =  —  from  (1.2),  (3.3),  (2.17)  we  have 


,  j  .2,  ,  T  n-2 

(grad  u)  dx  <  L  u 

|x-x0l<2 


(3.4) 
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Let  n  be  a  constant.  Multiply  (3,1)  by  the  function 

f c2 ( x) (u-n)  (u  >  n) 

\  o  (u  <  n) 


and  integrate,  we  obtain  similarly 

/  *;2{grad  u)2dx  <  /  (u-n)2(grad  C)2dx  •  (3.5) 

u>n  u>n 

We  denote  the  sets  u(x)  J  n  in  sphere  Jx-x^)  <  y  by  A  ,  B  respectively. 

^  n  $  u  n  f  m 


From  (3.2),  (3.5)  we  have 


f  (grad  u)2dx  <  L_  mes  A  — * — -  max 

A  .  5  n'“  ( 9 u ) 2  xeA 

n,u-tJp  r\,u 


max  [u(x)  -  nl 


Similarly  we  have 


(grad  u)2dx  <  L  mes  B  — - — r  max  (n-u(x)]- 

6  T1'W  (6w)2  XEB 


From  (3.4),  (3.6),  (3.7)  and  applying  the  result  of  [3J :  There  exist  constants  L,  and 
T  m  T  (Q,T)  (0  <  T ^  <  1)  such  that  for  any  x,  y  e  ft  we  have 

-T!  Ti 

|  u(  x )  -  u(y)|  «  L  r  |  x-y  |  .  (3.8 

/  xy 

Note.  In  (3]  the  inequality  (3.8)  is  proved  under  the  restriction  I x-y l  *  Lg rXy *  When 
| x-y 1  >  LgrXy,  the  validity  of  (3.8)  is  true  by  (2.17). 

We  now  turn  to  the  boundary  estimate.  From  (7)  we  see  that  there  is  a  local 


parametric  expression  for  a  little  part  of  f 


<2+V 


XX  =  XX(S . 5n-1)£C  ' 

in  the  neighborhood  of  any  point  e  T . 

Choose  variables  (£  )  suitably  such  that  inside  the  small  sphere  u>  with 

1  ,  _  .  n 

<2+V 

center  PQ  we  have  £  e  C  ,  on  point  pq  (^  , . . .  ,C^)  for  *  orthogonal  axis 

system,  and  £  coincides  with  the  direction  of  normal  on  u>  n  V .  The  method  for  choosing 

is  as  follows:  First  take  «  such  that  they  form  the  unit  orthogonal  axis 

1  n—  1 


system  on  PQ.  Penote 


cos ( N,  x  )  *  A. (£,...,£  )  (i  *= 

i  l  1  n-1 


on  T.  Then  the  inverse  function  of 

1  n 
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Ki  “  xi(C1 . en-1>  +  5n  i  V61  +  °15n . Vi  +  VlV*Y 


satisfies  all  our  requirements. 


From  (3.9)  we  have 


ajk  3x^  3x^ 


“jklr  “  Sic  (j'  k  *  1 . "> 

0 


i .  I  *0  (j=1,,..,n-1) 

^"'unr 


2  r  2  —  -  3*  . 

q  =  tu.  *  a  u  u  (u.  *  rr- 

i  ]k  3  k  3 


Equation  (1)  is  the  Euler  equation  of  the  variation  problem  6  /  F(q)dx  -  0  where 
F(q)  =  /  pqdq.  Under  the  transformation  we  have 

.  3  ( x  •  • . ,  x  )  . 

6  /  F,q)Jd5  -  o  fa  -  |___|) 


Jpa.  .  u ,  )  35  0 

3S  i3  3 

Differentiate  (3.12)  with  respect  to  £  we  have 

h 


3  3uh 

rr-UB..  +  B.  )  =  0 

13  3S  1 


B.  .  B  ptt.,  +  ^  tt,to  u,u 
1]  lj  q  ik  jm  k  m 


3  -  km  -  -  - 

B  .  *  T7~  (Ja  .  .  )pu  .  +  2JP  1  qa  .  .  tt — -  u  u  u 
i  3Ch  13  3  ij  3  *  m 


From  (3.10),  (3.14)  we  have 


8.1  *  (p<5  .  ,  ♦  —  u.u.)  • 

I3lp0  12  q  i  3  P0 

Hence  when  the  radius  of  u)  is  sufficiently  small/  form  (2.4),  (3.16)  we  have 


(3.17) 


I  ZXi  <  «  !  ZXl  • 

Denote  the  image  of  u>,  u)  n  ft,  w  n  T  by  w,  ft,  T.  (8)  becomes 

u  [  *  0 


(3.18) 


Let  C(S#  p 4  0)  be  the  function  defined  by  (3.3)  with  x,  x^  changing  to  £,  5q# 
restrict  in  this  case  that  the  sphere  with  center  £Q,  radius  p  lies  entirely  inside 


u>.  Multiply  (3.13)  by 


C2(u.  -  n)  (u,  >  n) 


(u,<  n) 
h 


and  integrate,  but  in  case  h  *  n  restrict  n  >  0.  Integrating  by  parts  we  have 

3u  3  2 

/  (J$.  .  57-  +  e.  )  TJ-  fC  (U.  -  n))dc  =  0  . 

~  r-  ^  \  13  3£ .  1  d€ .  h 

wnlUh^h)  3  i 


(3.19) 


Since  there  is  no  surface  integral  in  (3.19),  we  only  need  to  check  the  case  i  =  n. 
When  h  *  n,  on  f  by  (3.18)  we  have  u^  -  n  <  0.  When  h  ^  n,  j  *  n  from  (3.18)  we  have 

3\jl  .  3u  - 

=  0*  When  h  /  n,  j  ^  n,  from  (3.11),  (3.14)  we  have  3.1  =0.  And 

n  f  3Vf  n3!f 


the  b.  defined  by  (3.15)  has  the  property  b  )_  «*  0  when  h  ^  n,  the  reason  is,  from 
i  n  I  p 

*0  in  the  case  j  ?  n, 


3a 

(3.11)  we  have  a  I  •  0  (j  ^  n) ,  and  this  induces 

n3  I  f  dE> 


h  ?  n,  combine  with  (3.18)  we  have  b 


0. 


From  (3.19)  by  applying  Schwarz  inequality  and  (3.17)  we  have 


—  2  —  ^ 

i  (grad  u)  *  L_  mes(ft  n  A  ){ - -  max  [u 

ftnAv  ^  9  n#U  (0u)2  U«nAn  t  ^ 

k,p-9p  n,u 


[Zl  (U  -  M]  +  1}  .(3.20) 


Define 


VS . 


Vh . Vi'  -  V  (h  *  "> 

0  (  h  *  n) 

JO)  - 


(3.21 ) 


in  the  part  of  oi  outside  ft,  then  apply  (3.18)  we  have  u^  e  C  (~  )  W. ,(-,).  From 

(3.20),  (3.21)  we  have 
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f  (grad  u,  )2d£  4  2Lrt  mes  A  { — — -  max  [u  (£)  -  n]2  +  1} 

A  h  9  °'P  (0u)2  CeA  h 

n,u-0u  n ,u 


(restrict  n  >  0  when  h  -  n). 

By  the  same  process  we  can  obtain  two  formulas  similar  to  (3.4),  (3.7). 

From  the  result  of  [3]  we  have,  u^(h  *  1,...,n)  satisfy  a  Holder  condition  in  u> 

V 


i.e.  there  exist  L10  and  8)  (0  <  <  1)  such  that  for  any  x,  y  e  ft  n  u  we  have 


{ u ( x )  -  u(y) |  4  L1q |x-yj 


(3.22) 


}4.  Estimate  of  solution  in  the  neighborhood  of  ». 


Assume  that  the  solution  of  (1),  (8),  (9),  (11),  (2.17)  satisfies 

C(2)(H)  n  c(1)(fl  +  T)  . 

When  sP  £  E,  from  (2.5)  and  applying  the  Schauder  interior  estimate  twice  we  have 

•t  e  c(3+T,(fi)  n  c(1,<n  +  n  . 


(4.1) 


Let  0<8<1,0<8<1,X  is  a  positive  constant,  R,  R  are  large  positive 
numbers.  Let  r 


4(x)  -  C(x,  R,  5}  *  / 


I x |  R( 1+6 )  4  | X |  4  R(1-$  ) 

RX(1+0)X  R  4  |x|  4  R(  1+9  ) 


RX(1-0 .)X  Slisi  R(1-0  )  4  lx |  4  R  . 

’  V 


(4.2) 


Using  this  Q  we  can  obtain  (3.3)  also.  When  X  <  1  -  — ,  from  (2.17)  and  (4.2)  we  see 
that  (3.3)  is  also  true  when  C  *  C(x,  R,  «). 


Denote  the  inversion  point  of  x  about  the  unit  sphere  by  x,  then  (3.3)  becomes 
.2,„,4-2n,.  ,2_  „  .  r  2.  ,4-2n..  _2 


/  r  |xr^n(Axu)^dX  4  L2  /  u  lx|4"2n(AxC)  dx  . 

Take  X  »  2  -  n,  9  *  1,  R  *  •  and  let  ^  *  U ,  from  (4.2)  and  (4.3)  we  have 

/  (grad  u)2dX  4  L  j”  2 

|X|4*  X  11 


Similarly  we  have 


^  ( qrad  u)2dX  4  L  mes  A  - - — r  max  (u(X)  -  n)2 

A  f1  -1  *  2  °'M  (0U)2  XtA 

n,u(1+0)  ri  ,  p 


(4.3) 


(4.4) 


(4.5) 
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n  $  y(  1 +  9 ) 


2 

(grad^  u)  dX  <  mes  B 


r\,  y 


—  max  [u(X)  -  n 1‘ 


(9y)  XeB 

0  ,  U 


(4.6) 


where  A  ,  B  are  the  sets  of  u(X)  >  n  in  |x|  <  y  respectively* 
n ,  y  hf  u 

It  is  easy  to  prove  that  u(X)  e  in  a  neighborhood  of  the  origin  (which  includes 

the  origin).  From  (2.17),  (4.4),  (4.5),  (4.6)  and  applying  the  result  of  [3],  we  see  that 
constants  a  =  a(Q,  D  >  0  and  L14  exist  such  that  in  the  neighborhood  of  the  origin  the 
following  inequality  is  valid: 

i  at 


lu(X)  -  U1X=0I  <  l14'X' 


(4.7) 


where  u ! X=Q  =  u|  -  u(°°)  is  defined  from  (11).  Let  u(x)  -  u(«)  =  v(x),  then  (4.7) 


becomes 

I v ( x ) |  <  L14lx| . 

Apply  a  linear  transformation  x  =  x(y)  such  that  the  matrix  (a  ) 

i  j  x=® 

unit  matrix.  Then  (3.1)  becomes 

*1  (cii  Wj  '  0 

where  c. .(«)  -  5. . .  Combining  with  (4.8)  we  have 
13  13 


(4.8) 

becomes  an 

(4.9) 


lcij(y)  -  6i , 1  ‘  Li5lyl"a  • 


(4.10) 

Multiply  (4.9)  by  ^'v  and  integrate  to  obtain 

/  £2(grad  v)2dv  <;  j  v2(grad  £)2dy  (4.11) 

by  us'nq  the  same  method  as  (3.3)  was  obtained,  where  ^(y)  =  £(y,  R,  R)  is  determined 
from  (4.2)  by  changing  x  to  v.  Take  6  satisfying 

0  <  B  <  j  -  1,  0a*  (4.12) 

In  (4.2)  take  X  =  1  -  —  +  3,  R  =  °°,  9=1,  substituting  in  (4.11)  and  applying  (4.8)  we 
have 

f(r)  =  f  I y | 2  \  grad  v)2dy  <  L  r2  (  0  .  (4.13) 

I y I >r  17 

Estimate  bv  extending  the  method  of  [4].  Denote  the  unit  spherical  surface  with 
center  origin  by  ^ j,  let  r  <  R,  and  denote  lyl  =  y  we  have 
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f(r)  -  /  |y|2Xfgrad  v)Zdy  >  /  dw  /  (•— ) 


R  av>2  n- 1  +  2  X 


r< | v I 4R 


3u/ 


dp 


(4.14) 


>  f  [/  dp/  f  p  n+1  2Xdp]  du)  >  (n-2+2X)rn  2+2X  J  [v(  R)-v(  r )  ]2du 

a)  r  »■  a) 


where  in  the  above  expression  Schwarz's  inequality  was  used.  Let  R  -►  »  in  (4.14). 
Applying  (4.8),  (4.13)  we  have 

f  ( r )  >  (n-2  +  2X)rn"2  +  2A  f  v2(r)dui-  (n-2+2X)  /  lyl^VdS  .  (4.15) 


|y|-r 


Combining  (4.12),  (4.13),  (4.15)  we  have 

r 

ly  I  >r 

Let 


f  t  ,26-2  2  .  ^  _  2 ( 8-a) 

/  I Y I  vdy  <  L1Qr 


g  ( r )  *  /  1  y  I  2  Xc  .  .  -jp-  dy 

R  r <  |  y  1  4R  13  9yi  3yj 

Integrating  by  parts  and  applying  (4.9)  we  have 

2  X  3v 


Vr)  “  (/  -  /  )lyl  vcii  cos(N#  yi)ds  “ 

R  |yl-R  lyl-r  n  dyi  1 


(4.16) 


(4.17) 


(4.18) 


i  I  2 X— 1  3v  v  , 

ly  I  VC  T-—  cos  (N ,  y.)dy 

r<ly|ffi  13  3Yj  3 


where  N  is  the  normal  directed  toward  <».  When  R  *  °°  from  (4.10),  (4.13)  we  see  tht 
the  left  hand  side  of  (4.18)  has  a  limit,  from  (4.13),  (4.16)  applying  Schwarz  inequality 
we  prove  easily  that  the  last  term  in  the  right  hand  side  of  (4.18)  has  a  limit  when 
R  +  hence  when  R  *  » 

f  |y| 2Xvc  — —  cos ( N ,  y.)dS  (4.19) 

lyl-R  1  ^Yi  1 

has  a  limit.  If  this  limit  value  is  a  constant  p  ¥  0,  then  applying  (4.10),  (4.13), 


(4.15)  when  p  large  enough  we  have 

p2  4  L  (f  I y I  2 X f v (  {grad  v|ds)2  4  L  p  /  ly|2X  1  v2dS 

1R  I y ) =p  19  I y I =u 


f  | y | 2X( grad  v) 2dS 

I  y!  -u 


Integrating  we  qet 


4  (-L^pf  (p)f  *  (p)f 1  (p))/(n  -  2  +  2X) 


p2 ( n-2-2 \ ) 1np/L19 I R  4  -  ~  f2(p)|R  <  R> 


(4.20) 
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(4*13)  contradicts  (4.20)  when  R  +  «,  hence  the  limit  of  (4.19)  must  be  zero.  Let 


> 

i 


i 


* 


r  +  oo  in  (4.18)  and  denoting  g  (r)  by  g(r)  we  have 

,2X.  3v 


9<r)  *  -/  |yl^Avc  ~  cos(N ,  y  )dS  -  2 X  /  lyl2*”1  v  dy 


|y|=r 


ij  3y 


-2X  /  I  y  1 2  X  1  v(c  -  6.  .)  cos  (N ,  y  )dy  »  I  +  I  +  I 
)  y  J  >r  ij  i]  3yi  3  12  3 

Applying  (4.10),  (4.13),  (4.17)  we  have 


lyl  >r 


3N 


(4.21) 


g(r)  -  [1  +  0  ( r  a))f(r)  . 


(4.22) 


Applying  Schwarz  inequality  to  I1  and  applying  (4.10),  (4.13),  (4.14)  we  have 


I  <  r  J  lyl  v  dS  /  | y I  (c 

ly!=r  lyHr 


3v  2 

ij  W~  cos(N'  Yj^  ^  < 


nf  ( v) 
n-2+2X 


J  lyl  +  I V I  a0 ( | grad  v[)  dS  < 

,1  dN 

I  y  I  =r 


[1  +  0(r_a)]rf(r)/(n-2+2X)  J  |yi2X(grad  v)zdS 

I  y  I  — r 


-  (1  +  0(r  )) (rf(r)f'(r))/(n~2+2X) 


(4.23) 


Integrate  by  parts 


T  _  ,,  r  2 X+n-2  r  dv  ^  w  2\+n-2  ,  2 „  % 

12  ~2X  !  **  dy  /  v  —  dw  *  -X(y  /  v  d o>) 


3v 


2  X+n-2  f  2,  « 


+  X(2X+n-2)  /  y2X  n  J  v^du)  *  (1  -  ~  +  £)  f  |y( 
r  w  ly|*r 

v2dS  +  3(1  “  “  +  6)  /  lyl 2X  2  v2 dy  <  0 
I  y  I  -r 

the  validity  of  the  above  formula  has  used  (12)  and  the  relation 

,  2  X+n-2  r  2.  . 

( y  [v  dw)  =  0  . 

*  l_l=w 


2X-1 


From  (4.15)  we  have 


r  ,  ,2X-2-2a  2  ^  1 

J  I V I  v  dy  < 

lyl  <r 


n-2+2X 


ady  /  |y|  (grad  v)  dy 

ly  l*r 


<  L20r_2af(r> 


(4.24) 


(4.25) 


From  (4.10),  (4.13),  (4.25)  and  applying  Schwarz  inequality  to  estimate  we  have 

lTJ  *  L9i  /  ly|2A  1  “M  I  grad  v|dy  <  L  r -otf  (  r  )  .  (4.26) 

lyl >r  22 

Substituting  (4.22),  (4.23),  (4.24),  (4.26)  into  (4.21)  we  have 

[1  +  0  ( r”a)  1  f  (  r )  <  { r  1  +  0[r'a)l  (-rf(r)f'(r)),/(n-2+?\)}1/2 
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Integrating  we  get 


f(r)  <  [1  +  0(r  a)l (-rf' (r))/(n-2+2X) 


r  n-2+2A  P(r~a) 


f(r)  <  f<r0>(~)  e  0  <  L23r2"n"2X  (r  >  r0(C'  r)  > 


f  ,  ,2-n+2A  _  . 2  .  -28 

J  l Y I  (grad  v)  dy  <,  L  r 


I  y  f  >r 

Substituting  (4.27)  into  (4.15)  we  have 


Take  y  satisfying 


/  |y|_n+2Bv2dS  <  L  r" 

lyl-r  24 


0  <  Y  <  T  -  1,  Y<2B  . 

In  (4.2)  taking  X^l-j+y,  R«®,  0  «  1,  substituting  into  (4.1)  and  applying  (4.28) 


/  lyl2  n+2^(  grad  v)2dy  <  £2(grad  v)2dy  <  L  /  |y| 


-n+2y  2 


v  dy  <, 


_  f  -1-4S+2Y,  „  „  -4 8+2y 


Note  that  (4.29)  improves  (4.13). 

Repeating  the  above  process  we  prove  that 

f  lyl2  n+2^(grad  v)2dy  <  L  r  2y  . 

I  Y  I  >r 

Continuing  in  this  way,  after  a  finite  number  of  steps  we  have 

f  (grad  v)2dy  <  L  r2  n 

lyl>r  29 

We  need  the  following  result  of  [2]  extended  to  the  n  dimensional  case: 

f  ,  ,  ,2  _  ,  T  -2/n-1 

J  (grad  v)  dy  <  L  r  • 

I  y  I  >r  30 

Proof.  Let  S  be  a  closed  surface  which  lies  in  ft  and  contains  f  in  its 
interior.  Integrating  (1)  and  applying  (9)  we  have 


f  p  cos(N,  x.  )dS  »  f  o  ~  dS  «  n 

-  dX .  1  _  3N 


When  changes  to  x^  t  e,  the  changed  function  will  be  denoted  by  adding 


subscript  r.  Denote  bv  S  the  surface  5  translated  bv  (n,...,0,  t.  n  n) 

f. 

from  (4.32)  we  have 


/  (p  -r — )  cos(N,  x.)d§  «  f  p  cos(N,  x  )dS  *  0 

“  dx  E  1  dx  1  E 


S  i 


S  i 

e 


Differentiating  with  respect  to  the  parameter  e  and  let  e  »  0  we  have 


f  a  cos (N ,  x  )ds  *  /  aJ  cos(N,  x.)dS 

5  3xi  *  '5  W  3xt  J 


/  fp  ^  ^  ]  l  cos(N,  X  Ids  «  0. 

3e  axi^Ele«0  i 


After  linear  transformation  (S  changes  to  S)  we  have 

/  c  cos(N,  y  )dS  -  0  .  (4.33) 

S  ^  yi  ^ 

From  (4.30)  we  see  that  we  can  take  X  *  0  in  (4.21),  and  it  become  (by  using  (4.33)) 
q(r)  *  -  f  (v-v)c,.  cos  (N ,  y .  )  dS  (4.34) 

'  .  IT  .Kr  4  t 


v  *  J  vdS/mes  S 

lyl-r 


From  wirtinaer*s  inequality 


/  (v  -  V)  2cis  <  ~  I  t  ( grad  v)2  -  <|^)2ldS 

lYl-r  n'’  lyl-r  3N 

and  (4.10),  (4.34)  we  have 

/ri^T  [ 

9  r  *  2r  |y|,r(v-v)  dS  ♦  -g-  cos<N,  y^l^S  < 

!  f  ( qrad  v)2  -  (|^1  Ids 


+  L31r  a<qrad  V>21<1S  -  JTn^  U  +  0(r  a"  /  ( grad  v)2ds 

lyl-r  ■  ly  l=r 

=  -r  ( 1  +  0(r'c‘)]/(2/n^i)q'(r)  . 


Inteqratinq  we  have 


r0 , 2/n-1  0(r0  1 


g(  r )  <  g(  rQ  )  (~)  e 


(r  >  r  (p,  D) 


From  (4.10),  (4.17),  (4.35)  we  obtain  (4.31). 


Combining  (4.30),  (4.31)  we  obtain 


I  y  I  >r 

Substituting  into  (4.15)  we  have 


(grad  v)2dy  <  L^r*1  4  2  5  ( <$  >  0) 


lyl*r 


2  .  t  n-3-26 
v  ds  <  L33r 


Hence  we  have 


r  i  , 2-n+  6  2  r 

J  IyI  v  dy  -  J 


v  dy  -  J  r2  n+<^dr  / 

I  y  I  >R  R  I 

Reverting  to  the  variable  x  we  have 


y  l-r 


v2dS  <  L  R~^ 
34 


/  lyl 

|  x  |  >R 


2-n+S  2 .  -6 

v  dx  *  L35R 


(4.36) 


Let  where  U  is  defined  by  (11)  above,  then  from  (4.36)  we  have 

/  lyl  2  n+<^(grad  ^)  2dx  <  nL^R  ^ 


|x|>R 


Let  X  be  the  inversion  point  of  x  about  the  unit  sphere  and  let  -  *  U*  the  above 


formula  becomes 


I X 1  2  n  ^(cprad^  \Jj) 2dX  <  nL^u^ 


/  |X|-  "  "( qrad*,  \p)"dX  <  nL,cnu  .  (4.37) 

o  |x|<u 

Fixing  a  point  x  e  Q,  without  loss  of  generality  we  can  assume  temporarily  that 

-  0  x 


^(x°)  *  0,  hence  \p( x )  ■  «q  +  /  Q  v^dx^,  and  combining  with  (2.17)  we  have 


I  *(*)  I  <  LjglxJ 


I  *<X)  |  <  l36Ix 


-1 


(4.38) 


From  (4.37),  (4.38)  we  have 


e  w’  in  I X |  <  g0(Q,  D 


and  the  following  Sobolev  decomposition  formula  is  valid  [12]: 

I v| 1-nw  (x,  v)  |£-  dy 

I Y I <M0  I  Y 1 <UQ  3Yh 

where  £  and  w.  are  known  bounded  functions. 


<KX)  *  /  c<y)^(v)dy  +  / 


From  (4.37)  we  have 

f  I  Y|  1  n  |  grad  ^i|  dY  <  [/ 


(4.39) 


I  Y  |  <U 


I  v  |  <u 


|Y|6'ndY  J 


I  Y  | 


I  Y  |  <p 


2  n  ^(  grad  i|02dvj^2  . 


Hence  the  right  hand  side  of  (4.39)  ig  a  continuous  and  bounded  function  when  |x|  <  . 


Therefore  0(x)l  exists,  and 

x=» 


H>(  x  )  I  <  L 


3Q 


(4.40) 
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Denote  -  ^(«)  still  by  ^(x).  Then  from  (3.8},  (3.22},  (4.8},  (4.40}  we  obtain 


that  the  solution  of  (1),  (8),  (9),  (11  >,  (2.16)  is  bounded  in  E. 

All  the  requirements  of  Leray-Schauder  degree  theory  are  satisfied,  hence  theorem  2  is 
proved. 


§5.  Some  properties  of  subsonic  flow. 

Theorem  3.  The  solution  of  problem  (1),  (0),  (9),  (10)  is  unique  in  the  subsonic 
range. 

Proof.  From  (2)  we  obtain  that,  the  difference  v?  “  <^2^  -  ^  ^  of  two  solutions 


y?*2*  satisfies  the  following  equation 


.  (2)  _  (2 )  „  3  V  .  ,  & 

. . .  >  +  \  37;  - 0 

j  k 


where  u 


(1  ) 


(  1  ) 


(2)  _  3^ 


.  (u  ^  - 

ij  1  n  9x^3x^ 

(2) 


3x.  *  i 
i 


3x, 


and 


(5.1) 


b  =  [a  .  .  (  u 


(2) 


(2)  (1) 


-k  '"-,uk  '  Vi 


(1)  ,  .  ,  (2)  (2) 

...,u  )-a  (u  , 

n  ij  1  k  1  2(1) 

(1)  (1)  (2)  (1  ),  3  f' 

\  . v  n/tuk  ■  \  1  jtjt: 

i  3 


(5.2) 


From  (8)  we  have 


(5.3) 


From  the  last  section  we  have:  (after  subtraction  of  a  suitable  constant) 

satisfy  (2.6),  hence  we  have 

<?(")  -  0  .  (5.4) 

From  (5.1),  (5.2)  by  applying  the  maximum  principle  [5]  we  have,  when  is  not  a 

constant,  then  it  cannot  taken  a  positive  maximum  value  or  a  negative  minimum  value  in 

ft  •+■  T,  so  combining  with  (5.4)  we  have  yj  *  0.  This  proves  the  theorem. 

00  00 

Theorem  4.  There  exists  a  positive  constant  q^,  such  that  when  0  4  U  <  q^,  the 


I  1 


solution  of  problem  (1),  (8),  (9),  (10)  exists  in  the  subsonic  range  (this  solution  is 


unique  by  theorem  3).  Moreover,  the  function  Q  defined  by  (12)  satisfy 
Q(U)  e  C[0,  q")  ,  and 

lig  Q(U)  -  q  .  (5.5) 

U*qc-0 

The  proof  is  divided  into  several  points  as  follows: 

1)  Let  }  be  any  sequence  of  solutions  of  (1),  (8),  (9),  (11),  (12),  (13)  with 

n 

Q,  U  substituted  by  Qn,  Un«  If  sup  Qn  <  qc,  then  from  §2,  iii  we  have  that  liv^l  is 

uniformly  bounded  in  E,  {  ^}  is  compact  by  theorem  1,  hence  a  limit  function  exists, 

it  is  easy  to  prove  that  is  a  solution  of  (1),  (8),  (9),  (11),  (12),  (13)  also  and  its 

corresponding  Q  is  a  limit  point  of  the  sequence  {Q  } . 

n 

2)  Take  any  solution  \p  of  (1),  (8),  (9),  (11),  (12),  (13),  we  define  U  *  U(Q,yM 

by  (11).  Fixed  Q(0  <  Q  <  q^)  and  denote  the  infimum  of  all  U(Q,V>)  by  V(Q),  i.e.  a 

solution  sequence  }  exists  such  that  U(Q,v  )  *  V(Q).  By  1),  there  exists  a  limit 
n  n 

function  »  ^(Q),  which  is  the  solution  of  (1),  (8),  (9),  (11),  (12),  (13)  (where  we 
substitute  for  U  by  V  in  ( 1 1  ) )  . 

3)  Let  0  <  Qq  <  we  P^ove  that  the  interval  (0,  V(Qq)]  is  covered  by  the  set 

of  all  (U,  ifi)  (0  <  Q  <  Qq) . 

If  this  is  false,  i.e.  constant  Uq  exists  such  that  0  <  Uq  <  V(Qq)  and  the  problem 
(1),  (8),  (9),  (10)  (where  in  (10)  we  substitute  for  U  by  UQ ) ,  has  no  solution  in 
0  <  Q  <  Q0. 

If  there  is  a  positive  number  sequence  e  with  e  +  0(n  ♦  00 )  such  that  (1),  (8), 

n  n 

(9),  (10),  (where  in  (10)  we  substitute  for  U  by  Un  -  e  ) ,  has  a  sequence  of 

«  n 

solution  in  0  <  Q  <  Qq,  then  by  1),  the  limit  function  is  the  solution  of  (1), 

(8),  (9),  (10),  where  in  (10)  we  substitute  for  U  by  UQ,  in  0  <  Q  Q  ,  and  this 
contradicts  the  above  hypothesis. 

Hence  a  constant  U1  exists,  satisfying  0  <  <  Uq,  such  that  the  problem  (1), 

(8),  (9),  (10)  has  no  solution  in  0  <  Q  <;  Qq  when  U^  *  U  <  UQ. 

I  \ 


Now  apply  Leray-Schauder  degree  theory*  to  the  set  -»  <  U  <  V2  (UQ  +1^)  in  the 
space  E  (for  the  relation  of  the  element  ^  c  E  and  U  see  (2.6)).  We  obtain  a 
solution  to  problem  (1),  (8),  (9),  (11),  (12)  for  any  given  0  <  Q  <  and  this  solution 

possesses  the  following  property:  U  defined  from  (9),  (10)  satisfies  0  <  U  <  UQ*  Hence 

inf  U(Q0)  <  V(Q0)  . 

This  contradicts  the  definition  of  V(Qq),  which  proves  3. 

4)  The  function  V(Q)  is  strictly  monotone  increasing  in  the  interval  [0,  q^.). 

Otherwise  there  are  0  <  £>  <  Q_  <  q  such  that  V(Q  )  >  V(Q  ).  From  3)  we  have,  when 

12c  12 

0  <  Q  <  Q1  the  set  of  U(Q,v?)  covers  the  interval  [0,  VfQ^)],  hence  we  have 
Q0  e  [0#  and  a  such  that  U(Q0#  ^g)  “  V(Q2)f  contrary  to  theorem  3. 

V(o)  is  strictly  monotone  increasing  and  has  supremum  qc,  hence  V(qc  -  0)  exists. 

5)  Let 

<C =  v<<u  - °>  • 

c  c 

From  3)  we  have,  the  set  of  U(Q)  (0  <  Q  <  q  )  covers  the  interval  0  <  U  <  q°°.  In  other 

c  c 

words,  the  solution  of  problem  (1),  (8),  (9),  (10)  in  the  subsonic  range  exists  when 

0  <  U  <  q°°. 

c 

The  above  solution  is  also  unique  by  theorem  3,  hence  the  function 

Q(U),  0  <  U  <  q~  , 

is  determined  uniquely. 

For  fixed  U0  in  (0,  q*) ,  it  lies  in  the  number  set  U(Q(Uq),*£),  hence  we  have 

V(Q(Uq) )  <  UQ  .  (5.6) 


*  r 

From  §2,  iii,  but  substitute  Q  by  Qg ,  positive  constant  K  exists  such  that  the 
solution  of  (2.13)  satisfying  •  |  v  ||  <_K,  Take  F  be  the  part  of  sphere  j|tfj[  <_  K+ 1 
satisfying 

-w<  U  <  V2  (Un  ♦  Ut)  if  =  Uxt  +  f,,  D  ( « )  =  0) 

then  (2.13)  has  no  solution  in  the  neighborhood  of  the  boundary  of  F.  From  §2,  ii  we 
have,  the  learee  of  solution  is  1  in  F  when  *1  =  n,  combine  with  ^2,  i  and  applying 
Leray-Schauder  leqree  theory  we  obtain  that,  when  0  _<  ‘  ^  1 ,  (1)#  (8),  (9),  (11),  (12) 
in  F  has  at  least  one  solution. 


1 


6) 

From  the  definition 

of  q°°  there  exists 
c 

% 

+  q  -  0,  such 
c 

that 

V<2  >  ♦  q” 
n  c 

If 

U  satisfies  0  <  u 

<  V(0  )  and  Q(U) 

'  n 

> 

V 

Combining  with 

(5.6) 

we  have 

u  >  v(Qn: 

)  ,  this  contradicts 

U  <  V(Q  ).  Hence 
n 

when 

0  <  U  <  V(Q  ). 

n 

we  have 

0  <  Q(U)  <  Q  <  q  .  (5.7) 

n  c 

7)  Let  {U  }  satisfy  0  <  U  <  V(Q  ) ,  U  *  U  ,  Q(U  )  ♦  Q  .  From  (5.7)  by  applying 
m  m  n  m  o  m  u 

1)  we  obtain  that  the  limit  function  ^  is  the  solution  of  (1),  (8),  (9),  (10),  where  U, 
Q  are  given  by  UQ,  Q0  in  (10),  (12)  respectively.  Applying  theorem  3  we  have 
Q0  *  Q(Uq),  in  other  words  Q(U)  e  C[0,  V(Q^)].  Letting  n  ♦  »  we  have 

Q(U)  e  C[0,  q“)  . 


8)  From  the  above  points  we  see  easily  that  lim  Q(U)  »  q  .  then  the  functional 

U4q°°-0  C 

c 

value  of  the  continuous  function  Q(U)  oscillates  finitely  in  the  interval  (Qo'qc*'  an<3 

every  point  in  (Qo'qc^  *s  the  P°*nt  of  Q(U)  when  U  -►  qc  -  0.  Hence  we  have 

U  ♦  q  ,  U*  ♦  q  such  that  Q(U  )  ♦  ,  Q(U*  )  Vo  (Q„  +  q  ).  By  1),  the  limit  functions 

mcmc  mum  uc 

^ ,  V  of  corresponding  solution  sequences  },  &  * }  are  solutions  of  (1),  (8),  (9), 

m  m 

(11),  (12),  (13),  where  Q  is  given  by  Qg  and  V2  (Qg  +  qc>  respectively.  Hence 

-  V«*o 

U(Q0)  -  qc  *  U( — 5 — ) 

00 

or  the  problem  (1),  (8),  (9),  (10)  (substitute  U  by  )  has  at  least  two  solutions  in 
the  subsonic  range,  which  contradicts  theorem  3,  hence  (5.5)  is  true. 

The  theorem  is  proved  completely. 

Condition  (11)  restricts  the  incoming  flow  is  in  the  positive  direction.  Now 

remove  this  restriction  by  only  assuming 

oo  00  00 

u  *  (u. , . . . ,U  )  (5.8 ) 

i  n 

is  given.  Because  (1),  (8)  is  invariant  under  axis  rotation,  hence  by  theorems  3,  4  we 
have,  the  solution  of  (1),  (8),  (9),  (5.8)  exists  and  is  unique.  Let 


*  =  ^ 

and  regard  ^  as  a  solution  of  the  linear  nroblem,  by  theorem  1  we  have 


(5.9) 


H  i|ill 


n-  2 , 2  +  t 


39' 


(5.10) 
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Now  prove  that  the  solution  of  (1),  (8),  (9),  (5.8)  depends  continuously  on  u  by  the 
following  sense* 

Theorem  5.  There  exists  a  positive  constant  L40'  such  that  for  anv  two  vectors 
00  2 

u  '  in  the  class 


,  00,  O  0*  Q  CO 

U  <  Vi-*—  U,,...,-* -  U  ) 

lu  I  |u  I 

(the  definition  of  function  V  see  theorem  4)  the  corresponding  solutions  of 

(1),  (8),  { 9 )  ,  (5.8)  satisfy 


(2) 


.(1) 


‘  L«|u”'2  -  u”'1'  (5*n) 
where  the  relation  of  y  and  4)  is  given  by  (5.9). 

Proof.  Let  $  *  i|/2*  -  4/ 1  *  ,  from  (2)  we  obtain  that  £  satisfies  (5.1),  (5.2). 

Prom  (8)  we  have 


,  00, 2  •»,  1 

wlr  =  ~(ui  ”  ui  ’  cos(N'  Vr  • 

Applying  (6),  (5.10)  it  is  easy  to  show  that  the  b^  defined  from  (5.2)  satisfy  conditions 
(1.7),  (1.55),  (1.59),  and  from  (5.10)  we  get  (5.4).  Hence  from  theorem  1  we  have  (5.11). 
This  completes  the  proof. 

Theorem  6.  For  every  non  zero  subsonic  flow  around  a  given  profile,  cannot  be  taken 

on  >1  or  at  ».  In  other  words,  q^x  can  only  be  taken  on  T.  And 

<W<L41|U">  *  (5-12> 

Proof.  If  q  >0  is  taken  by  p  c  Q,  then  after  a  rotation  of  axis  we  arrive  at 

grad  u|p  -  (qmaxr  0,...,0)  and  (1),  (8)  remain  unchanged,  hence  u^  takes  the  maximum 

value  at  P.  And  u1  satisfies  the  elliptic  equation  (3.1),  hence  from  the  maximum 

principle  (51  we  have,  is  a  constant,  i.e.  u1  =  amax,  hence  u^  -.*.=  un  *  0*  Prom 

(8)  we  have  n  „  =0*  This  is  a  contradiction, 

max 

Prom  (5.10)  and  applying  the  result  of  (2],  we  have  that  qmax  cannot  be  taken  at 


CO  2 

(5.12)  is  a  special  case  of  (5.11)  with  u  =  0.  The  proof  of  the  theorem  is  thus 
complete . 

Theorem  7.  u)  exists  for  any  subsonic  flow  in  the  neiqhborhood  of  »  *nd  (4.8)  is 

valid.  An  the  subsonic  flow  in  the  whole  space  is  uniform* 


-32- 


Proof.  For  subsonic  flow  in  the  neighborhood  of  •,  because  (4.4),  (4.5)  and  (4.6) 

are  valid,  hence  the  existence  of  u|  and  the  validity  of  (4.8)  follow  from  [3). 

x*® 

Concerning  subsonic  flow  in  the  whole  space,  u|  ^  exists  by  the  above  argument, 
hence  it  is  a  special  case  of  flow  around  profile,  i.e.  (4.32)  is  true  for  every  closed 
surface.  Hence  (5,t0)  is  valid,  and  we  can  apply  theorem  6,  i.e.  except  q  is  a  constant, 
q#  cannot  be  taken  by  any  finite  point  and  ®.  And  when  q  is  a  constant,  u  is  a 
constant  vector  by  theorem  6,  i.e.  the  flow  is  uniform. 

The  theorem  is  thus  proved. 
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